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ABSTRACT 

We propose a framework for separating the information contained in the CMB 
multipoles, ae, n , into its algebraically independent components. Thus we cleanly 
separate information pertaining to the power spectrum, non-Gaussianity and pre- 
ferred axis effects. The formalism buil ds upon the r ecentl y proposed multipol e vec- 
tors (iConi. Huterer fc Starkmanll^ ISchwarz Fall^nHKatz fc WeeksiElfil. and 
we elucidate a few features regarding these vectors, namely their lack of statistical 
independence for a Gaussian random process. In a few cases we explicitly relate our 
proposed invariants to components of the n-point correlation function (power spec- 
trum, bispectrum). We find the invariants' distributions using a mixture of analytical 
and numerical methods. We also evaluate them for the co-added WMAP first year 
map. 
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1 INTRODUCTION 

The remarkable quality of the WMAP data iBennett fc all 
2003) has led us to a new era in observational cos- 
mology. Yet, various claims for "unexpected" non- 
Gaussia n signals cast a shadow o n the p urity of the 
data JCoDi^hiterer^k Starkman 2003; Erikseii^^i 

f; boles fc all 120031: ICruz fc all 120031: Hansen fc a 
lHansen. Bandav fc Gorskil Eool IVielva fc a 
iKomatsu. Spergel fc Wandehl 120031: IParkl \ 200A 
Larson fc Wandeld 12004 iMukheriee fc Wand |2004 
Eriksen fc alll2004al : ISchwarz fc al2004 : IOliveira-Costa fc all 
2004lh ~Given that there is no plausible theoretical explana- 
tion for these signals, a natural worry is that the maps are 
not, after all, fully free from systematics errors or galactic 
contamination. 

Most notably several groups have reported evidence for 
a preferred axis being selected by large-a ngle fluctuations, ei- 
ther i n the form of multipole planarity (lOl^eira^osta^^d 
| 2004l: ICopi. Huterer fc Starkman! 120031: ISchwarz fc alll200l 
Cole s fc a ll l2003ft . North-South asymmetries in the power 
spectrum, three-point function, or bispec t rum ([Eriksen fc~a 
2004al: lHansen. Bandav fc Gorskil 12004 lLand fc Maeueiic 
2004 ;), as well as u sing other methods feriksen fc ajf l2004b: 
Hansen fc all2004h . 

In assessing these asymmetries it is important to dis- 
tinguish issues of non-Gaussianity (which should be rota- 
tionally invariant) from those of anisotropy (existence of a 
preferred axis). Apart from a subtlety in the definition of 
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statistical ensemble jFerreira fc Magueiidll997h these issues 
should be clearly separated. Unfortunately no systematic ap- 
proach for extracting all the independent invariants under 
rotations from a given se t of multipoles, ai m , has been pro - 
posed. The formalism of ICopi. Huterer fc Starkman! (|2003T> 
does not prod uce invariants. The invariant n-point correla- 
tion function jMagueijdll995l:l Ferreira^Maguciio fc Gorskil 
I199&1 iMaeueiiol l2000l: iMaeueiio fc Medeirosl l2003ft . on the 
other hand, is awkward to apply and often contains redun- 
dant information. 

In this letter we remedy this deficiency in the current 
formalism. 



2 STATEMENT OF THE PROBLEM 

Multipoles are irreducible representations of SO (3), so their 
21 + 1 degrees of freedom should split into 21 — 2 invariants 
and 3 rotational degrees of freedom. Ideally we would like 
to process the {a^ m } for a given mutipole I into the power 
spectrum Ci (the Gaussian degree of freedom), 2^ — 3 in- 
variant measures of non-Gaussianity, and a system of axes 
(assessing isotropy). In addition one should build 3 invari- 
ants per multipole pair, for example the Euler angles relating 
the two sets of multipole axes. The latter encode inter-scale 
correlations. 

No o ne has eve r acco mplished this project for general 
£, but see lMagueiiol Bl995i) for a Quadrupole solution. Mul- 
tipoles are equivalent to symmetric traceless tensors of rank 
£. The problem is then to extract the independent invari- 
ant contractions of these tensors plus a system of axes - 
a generalisation of the concept of eigenvalues and eigen- 
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vectors. Typically the invariants produced by this formal- 
ism are related to the n— point correlation function (bispec- 
trum, trispectrum, etc.) The formalism becomes very com- 
plicated very quickly, and no systematic breakdown of inde- 
pendent n— point correlation function components has ever 
been achieved. 

Alternatively one may extract from the {ai m } a length 
scale and £ inde pendent unit vectors (the multipole vectors) 
as proposed bv | S chwarz fe"all <l2004l) , and discussed further 
bv lKatz fc Week! l20o3) ~This approach is considerably sim- 
pler and unsurprisingly it has been taken further. However 
we see a (correctable) problem with this approach. The mul- 
tipole vectors are not rotationally invariant, and so they mix 
up the issues of isotropy and non-Gaussinity. We propose to 
correct this shortcoming by taking an appropriate number 
of independent inner products between these vectors. These 
are the sought-after invariants, and in some simple cases we 
relate them to the n— point correlation function. We will 
also extract from the multipole vectors a system of axes, en- 
coding the multipole directional information. It is these axis 
that are to be used when testing isotropy. 



3 THE EXAMPLE OF THE QUADRUPOLE 

We start by considering the quadrupole, which as shown 



We start by 
in lMagueiiol 

8T 2 = Q i jX i x i 



1995) may be written as 



(1) 



where Qij is a symmetric traceless matrix and x 1 are carte- 
sian coordinates on the unit sphere. From this matrix one 
may extract three eigenvectors and two independent combi- 
nations of invariant eigenvalues Xi. These are essentially the 
power spectrum C2 (related to the sum of the squares of Xi) 
and the bispectrum B222 (related to the determinant of the 
matrix, or the product of Xi). 

In contrast, following the formalism of iKatz fc Week^ 
( 2004) one writes 



STo = Ao 



Li Lj — —SijL 



(2) 



where A is a scale and {L],L 2 } are two units vectors, en- 
coding the information on non-Gaussianity and anisotropy. 

The vectors {Lj,L 2 } are not invariant, but one can con- 
struct a n invariant by taking X = L 1 ■ L 2 . As pointed out 
to us bv lStarkman fc Schwarz) j2004l) . one may easily check 
that the eigenvectors of Q are 



V 1 = L 1 + L 2 
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These have corresponding eigenvalues A(X + 3)/ 6, A(X — 
3)/6, and —AX/ 3, respectively. Using results in iMagueiiol 
(1995) one may therefore prove that the power spectrum 
and bispectrum are: 
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Figure 1. The average distribution of the modulus of the 2£ — 3 
dot products Xij (with the two anchor vectors chosen at random) 
for multipoles 1 = 2 — 10. The dotted line corresponds to the ana- 
lytical expression found for the quadrupole. We have also plotted 
the invariant values for the WMAP first year data, with the an- 
chor vectors chosen at random (short lines along the bottom of 
each panel). 
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These formulae bridg e the two formali sms. 

It was proved in IMagueiiol (|l995h that for a Gaussian 
process C2, I2 and the eigenvectors are independent random 
variables, and that I2 is uniformly distributed in the range 
[—1,1] (and C2 has a xi distribution). We thus can prove 
that X is statistically independent from C2 (but not from 
A) and from the eigenvectors. By directly evaluating the 
Jacobian of the transformation we find that its distribution 



P(X) = 27 



1 — X 



(X 2 +3) 5 / 2 



(9) 



so that the normalised bispectrum is 



that is, it is not uniform. We have confirmed this result with 
Monte-Carlo simulations (see Fig. top left panel) . 

This elucidates an interesting feature of multipole vec- 
tors. Even though they are algebraically independent (in 
the sense that they contain no redundant information given 
a concrete quadrupole realization) they are not statisti- 
cally independent. If vectors L 1 and L 2 were statistically 
isotropic (which they are) and statistically independent, 
then X would be uniformly distributed. As © shows this is 
not the case; hence vectors L 1 and L 2 are statistically cor- 
related. Specifically they prefer being orthogonal to being 
aligned. 



Multipole invariants and non-Gaussianity 3 



: i 




i ^ 

i h - 


-, , , i 1 


i 1 

fHJ 



0.2 0.4 0.6 0.8 1 

t 6 




0.2 0.4 0.6 0.8 

i=B 



0.2 0.4 0.6 0.8 1 

1=9 




0.2 0.4 0.6 



0.2 0.4 0.6 0.8 1 




Figure 2. The distribution of the modulus of the dot products 
for multipoles 1 = 2 — 10, with anchor vectors L 1 and L 2 chosen 
as the most orthogonal pair of multipole vectors. The solid lines 
are for dot products Xu and X2i, and the dotted line shows the 
X12 distribution. The WMAP results are also plotted (short lines 
along the bottom, dashed line the X±2 result). 



4 THE GENERAL PROCEDURE 

For a general multipole we have £ vectors and we could take 
as invariants all possible inner products between them. This 
would clearly lead to much redundant information, in con- 
tradiction with the requirements laid down in Section |5] A 
possible way out is to select two "anchor" vectors L 1 and 
L 2 , consider their dot product X12 — L 1 ■ L 2 , and then for 
3 ^ i ^ £ the two products Xu = L 1 ■ U and X 2 i = L 2 ■ V . 
We thus obtain 2£ — J> algebraically independent invariants. 
In Fig. 0we plot their distribution for 2 ^ I ^ 10 after the 
order and direction (±) of the £ vectors has been randomised. 
The anchor vectors L 1 and L 2 are therefore selected at ran- 
dom, so that all Xij invariants for a given multipole have 
the same distribution. We checked that the invariants X are 
uncorrelated. As we can see this distribution is £ dependent, 
and the tendency for multipole vectors to seek orthogonal 
directions is less pronounced for higher multipoles. We have 
also plotted these invariants as computed for the WMAP co- 
added masked first year map (details in lLand fc Maeueiicl 
(2004)). No evidence for non-Gaussianity is found, but we 
defer a closer scrutiny to a future publication. 

The proposed procedure provides an interesting Gaus- 
sianity test. In particular it can be easily applied to extract 
all the relevant information for high £; in contrast isolat- 
ing such information from the n-point correlation function 
is most cumbersome. However the suggested algorithm suf- 
fers from the drawback that the anchor vectors L 1 and L 2 
are selected at random, and so the procedure is not repro- 
ducible. 

We may correct this by imposing a criterium for their 



selection, for example L 1 and L 2 may be taken as the most 
orthogonal among the given £ vectors. We may then form 
the "eigenvectors" (three orthogonal vectors): 

V 1 = L 1 + L 2 (10) 
V 2 = L 1 — L 2 (11) 
V 3 = L 1 x L 2 (12) 

as the natural variables for encoding the multipole direction- 
ality. For invariants we may take the power spectrum and 
the 2£ — 3 quantities: 

L 1 ■ L 2 (13) 
L 1 ■ V (14) 
L 2 ■ V (15) 



X 12 

Xu 

X2i 



(for 3 ^ i ^ £). These encode all the relevant non-Gaussian 
degrees of freedom. 

The transformation from atm into these variables is 
invertible (up to discrete uncertainties related to branch 
choice) and provides a solution to our problem, as phrased 
in Section [5] The proposed variables for a given multipole 
£ are the power spectrum Ce (Gaussian degree of freedom), 
the 2^ — 3 inner products X (non-redundant non-Gaussian 
invariants), and the orthogonal vectors {V 1 , V 2 ,V 3 } (mea- 
sures of anisotropy) . The procedure reduces to the one found 
for the quadrupole when £ = 2. Within this framework the 
inter-£ correlations are measured by the Euler angles relat- 
ing the systems of axes {V 1 ,V 2 , V 3 } associated with each 
pair of multipoles. These should be uniformly distributed for 
a Gaussian distribution or indeed for any theory in which 
the various £ are uncorrelated. 

In Fig. |2]we plot distributions for the X invariants with 
anchor vectors L 1 and L 2 defined as the most orthogonal. 
We have used 12,500 realizations to make these histograms, 
f (^12) peaks around zero. The other X distributions are 
the same. We also plotted the invariants for the WMAP first 
year data, again finding no evidence for non-Gaussianity. 
We have checked that the eigenvectors, V™, are uniformly 
distributed. 

Notice that the invariants X cannot be independent 
variables, since their ranges of variation are interconnected. 
This is to be compared with the lack of independence among 
the various vectors U , as demonstrated in the previous Sec- 
tion. 

Naturally we could have defined the two anchor vec- 
tors L 1 and L 2 in different ways, for example the two most 
aligned vectors. The invariant X12 would then be peaked 
around one. We plot the counterpart of Fig. [5] with this al- 
ternative definition in Fig. 



5 BRIDGING THE TWO FORMALISMS 

For higher multipoles, relating the proposed formalism and 
the n— point correlation function, as done in Section [!|] for 
the quadrupole, becomes very involved. For instance, for the 
octopole the counterpart of Q and is: 



ijk 
XX 



(16) 



8T 3 = Qijkx % x J x k = A 3 (^LIL 2 lI - ^SijRk^jx 

where one should note that terms in {i,j, k} are sym- 
metrized. The remainder Rk is: 



4 Kate Land and Joao Magueijo 




- n 


H || i | | i 

i. - 
J j 

i : 

: 


: 


1_ i- 

i ih-L , ,- 



0.2 0.4 0.6 0.8 1 

1=7 




0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 

1=9 



0.2 0.4 0.6 0.8 1 

(=10 




0.2 0.4 0.6 



Figure 3. The distribution of the modulus of dot products for 
multipoles £ = 2 — 10, with anchor vectors L 1 and L 2 chosen 
as the most aligned pair of multipole vectors. The solid lines are 
for dot products Xn and X%i, and the dotted line shows the 
X\2 distribution. The WMAP results are also plotted (short lines 
along the bottom, dashed line the X\q. result). 



R k = h {1 ■ L 2 Lf 



(17) 



components {ai m } can be simply separated into a set of 
algebraically independent invariants - the power spectrum 
and a set of 21 — 3 non-Gaussian statistics - and a set of 
axes encoding the multipole directionality. The Euler angles 
relating sets of axes associated with pairs of multipoles mea- 
sure inter-^ correlations. This provides an elegant answer to 
a long unsolved problem - how to process the information 
contained in a given map into its relevant non-redundant 
degrees of freedom. 

Even though we have computed some of the proposed 
invariants for the WMAP first year map, we defer to a future 
publication a more systematic application of this framework 
to real data. We stress that the formalism is easily appli- 
cable for high I multipoles. This is to be contrasted with 
the n-point correlation function. E ven though it is trivial to 
evaluate the bispectrum at high I l|Land fc Magueiioll2004l) 
it becomes very difficult to distill all the non-redundant in- 
formation contained in a given multipole in terms of suffi- 
ciently high order components of the correlation function. 
The proposed formalism is far more efficient. 
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Again, Tr{Q 2 ) = C 3 , and therefore one can show that: 



C 3 



where 



6-9 



(9 + 7S 2 + 6S 2 ) 



Si 

S 2 



(L 1 ■ L 2 ) 2 + (L 2 ■ L 3 ) 2 + (L 3 ■ L 1 ) 2 
(L 1 ■ L 2 ) x (L 2 ■ L 3 ) x (L 3 ■ L 1 ) 



(18) 

(19) 
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For higher multipoles 

Q<i...<< = A t (Ll...L% - |% 1!2 (L (1 -i 2 4)-^) ( 21 ) 

The components of the n-point correlation function may 
be obtained from the various contractions of Q. They are 
clearly a function of the dot product of the various multipole 
vectors. However, the formulae become progressively more 
complex to derive. With the introduction of our alternative 
procedure we suggest that these complicated formulae are 
not necessary. The X variables provide a better basis for 
describing the data. 



6 CONCLUSIONS 

Harmonic (Fourier) space is the natural arena for comparing 
theory and observation for Gaussian theories. In several past 
studies it was also found to be a useful ground for testing 
the hypothesis of Gaussianity. In this paper we showed fur- 
ther that the degrees of freedom in the spherical harmonic 
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